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I . INTRODUCTION 


Aerodynamic  or  inertial  asymmetry  of  a fin-stabilized  non- 
rolling projectile  results  in  the  round's  trimming  to  a fixed 
attitude  once  launch  disturbances  damp  out.  This  preferential  trim 
produces  a continuously  increasing  divergence  of  the  flight  path  from 
the  intended  trajectory.  To  avoid  this  effect,  the  projectile  is 
given  some  small  roll  rate,  either  in-bore  or  in  free  flight.  The 
trajectory  of  a rolling  projectile  still  suffers  an  aerodynamic 
jump  due  to  the  asymmetry,  but  the  continuous  divergence  is  eliminated. 
Murphy  and  Bradley1 2  have  demonstrated  that  the  magnitude  and  direction 
of  the  aerodynamic  jump  are  related  not  only  to  the  steady-state  roll 
rate  but  also  to  the  variation  in  roll  rate  between  launch  and  steady- 
state  values.  Their  analysis  assumes  that  the  projectile  enters  free 
flight  with  roll  rate  equal  to  its  in-bore  value.  This  assumption 
neglects  the  perturbation  of  the  roll  due  to  transit  of  the  muzzle  blast. 

To  produce  a rolling  moment  on  a fin-stabilized  projectile, 

differential  fin  cant  is  generally  used.  In  the  reverse  flow  of 

the  muzzle  blast,  the  direction  of  the  rolling  moment  is  opposite  to 

that  in  free  flight.  As  a result,  the  roll  rate  decreases  from  the 

launch  value.  For  projectiles  fired  from  smoothbore  guns,  a reverse 

2 

spin  may  be  produced.  Fansler  and  Schmidt  extended  the  computations 
of  Murphy  and  Bradley1  to  include  the  case  of  projectiles  spinning  up 
to  a positive  free  flight  value  from  an  initially  negative  roll  rate. 
Additionally,  Fansler  and  Schmidt  calculated  the  transverse  linear  and 
angular  momentum  imparted  by  the  muzzle  blast  to  a projectile  with  an 


1.  C.  H.  Murphy  and  J.  W.  Bradley , "Jump  Due  to  Aerodynamic 
Asymmetry  of  a Missile  with  Varying  Roll  Rate3  " BRL  R 1077 3 
U.  S.  Army  Ballistic  Research  Laboratories  3 Aberdeen  Proving 
Ground^  Maryland,  May  1959 1 AD  219312. 

2.  K.  S.  Fansler  and  E.  M.  Schmidts  " The  Influence  of  Muzzle  Gas- 
dynamics  upon  the  Trajectory  of  Fin-Stabilized  Projectiles, , " 

BRL  R 1793,  U.  S,  Army  Ballistic  Research  LaboratorieSs  Aberdeen 
Proving  Grounds  Marylands  June  1975s  AD  B005379L. 

7 


aerodynamic  asymmetry.  In  a sample  case,  they  demonstrated  that  the 
increase  in  jump  due  to  negative  roll  rate  is  cancelled  by  the  component 
of  jump  due  to  the  projectile  transverse  velocities  imparted  by  the 
muzzle  blast.  However,  since  the  integration  of  the  equations  of  motion 
was  performed  numerically,  a closed-form  solution  was  not  obtained,  and 
no  general  statement  regarding  this  cancellation  effect  could  be  made. 

The  present  report  used  an  analytical  approach  to  obtain  a closed- 

form  expression  for  the  jump  of  a slightly  asymmetric,  fin-stabilized 

projectile  due  to  muzzle  blast.  The  flow  is  modelled  using  a quasi- 

2 

steady  approximation  which  provides  for  direct  evaluation  of  the  gas- 

dynamic  forces  on  the  projectile  during  launch.  From  the  resulting 

loading  history,  the  motion  of  the  projectile  through  the  muzzle  blast 

and  into  free  flight  is  determined.  Once  the  initial  dynamics  of  the 

round  are  known,  its  subsequent  trajectory  may  be  obtained  by  numerically 

1 2 

integrating  the  equations  of  motion  * . In  the  present  report,  the 
expression  for  the  trajectory  deflection  due  to  asymmetry  is  expanded 
as  a Taylor  series  and  compared  to  the  results  of  the  numerical 
integration.  Substitution  of  this  series  into  the  jump  equation  pro- 
vides a general  prediction  of  the  muzzle-blast-induced  jump  of  asymmetric 
projectiles. 


II.  ANALYSIS  OF  GAS DYNAMIC  LOADINGS 

At  launch,  a fin-stabilized  projectile  is  subject  to  gasdynamic 

loadings  in  two  distinct  muzzle  flows:  the  in-bore  expansion  and 

the  muzzle  blast.  Figure  1.  In  this  report,  the  analysis  of  Fansler 
2 

and  Schmidt  is  used  to  describe  the  muzzle  flow.  They  modelled  the 
in-bore  flow  as  a one-dimensional,  unsteady  expansion  and  the  external 

3 

flow  as  a quasi-steady  , axially  symmetric  jet  flow.  Transverse  gas- 

3.  K.  Oswatitsch , " Intermediate  Ballistics , " Deutsche  Luft  und 
Raumfahrt  FB  64-37 , DVL  Bericht  358 , December  1964 3 AD  473249. 
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Figure  1.  Sources  of  Gasdynamic  Loadings  During  Launch 
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4 

dynamic  loads  are  assumed  to  be  generated  only  upon  fin  surfaces  . 

These  loadings  are  computed  from  the  local  flow  properties  and  two- 
dimensional,  thin  airfoil  theory  without  correcting  for  flow  inclination, 
tip  effects,  or  wing-body  interference.  These  approximations  eliminate 
the  influence  of  fin  and  projectile  geometry  and  produce  an  upper 
bound  upon  muzzle  blast  loadings. 

Using  this  model,  the  loadings  upon  the  fins  may  be  integrated 
through  the  muzzle  blast  to  obtain  an  expression  for  the  momentum 
transferred  to  the  projectile.  Fansler  and  Schmidt  determined  non- 
dimensional  impulse  functions  that  are  independent  of  the  projectile 
and  flow  geometry  and  vary  only  with  propellant  gas  Mach  number  prior 
to  shot  ejection,  Vp/c^,  Figure  2.  Here  is  the  projectile 


Figure  2.  Comparison  of  In-bore  and  Muzzle-Jet  Impulse 
[for  d/D  = 1] 


4.  W.  Gretler,  Intermediate  Ballistics  Investigations  of  Wing 
Stabilized  Projectiles , " German  Air  and  Space  Research  Report 
67-92 , FSTC-HT- 23-22-69- 72 1 1967. 
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velocity  and  Cj  is  the  prelaunch  propellant  gas  speed  of  sound.  The 

total  impulse  is  the  sum  of  the  in-bore  and  muzzle  jet  function  values, 

P * P + (d/D)  P,  . Here  d is  the  distance  from  obturator  to  fin  and 
o b 

D is  the  gunbore  diameter.  This  function  is  used  to  calculate  the 
transverse  angular  and  linear  velocities  and  reverse  spin  imparted  to 
the  projectile  during  transit  of  the  muzzle  gases. 

The  coordinate  systems  used  in  this  report  are  standard  in 
ballistic  calculations1,  Figure  3.  The  nonrolling  coordinates  are 

Xg ,Ye,Z(  - EARTH  FIXED  SYSTEM 

X,Y,Z  - MISSILE  FIXED, NON-ROLLING  SYSTEM 

3,  $ - ANGLES  OF  ATTACK  AND  SIDESLIP  RESPECTIVELY 


a 


Figure  3.  Coordinate  Systems 

used  to  describe  the  projectile  angular  motion,  while  the  earth-fixed 

coordinates  are  employed  in  the  solution  for  projectile  jump.  At 

separation  from  the  gun,  the  projectile  velocity  vector  is  assumed  to 

lie  along  the  X axis.  Since  it  is  of  interest  to  examine  muzzle 
0 

blast  effects  caused  by  asymmetry,  we  assume 
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where  the  subscript  m referes  to  values  immediately  at  separation. 

Here  the  prime  superscript  denotes  rate  of  change  of  the  variable 
with  respect  to  distance  in  calibers  traveled  by  the  projectile. 
Further,  we  restrict  the  analysis  to  projectiles  launched  from  smooth- 
bore guns: 

/ = 0 . 
m 

For  simplicity,  the  aerodynamic  asymmetry  to  be  treated  is 
postulated  as  two  opposing  fins,  inclined  at  an  angle  e (in  the 
same  sense  for  both  fins)  with  respect  to  their  normal  orientation. 

The  differential  fin  cant  angle,  2 6,  Figure  4,  between  opposing  fins 


AXIS  TO  FIN  C.P. 


Figure  4. 


Illustrating  Differential  Fin  Cant 


is  unaltered  and  the  steady  state  roll  rate,  is  maintained.  The 

transverse  linear  velocity  imparted  to  the  projectile  in  the  muzzle 
2 

flow  is  computed  to  be: 


# * i r = (T  * 1)  p 2A,  £ J P (e  e (1) 

P P Mp  Vp 

where 

<}>e  = initial  orientation  of  asymmetric  force  in  free  flight. 

vq  = muzzle-blast  imparted  velocity  in  Yg  direction 
wq  = muzzle-blast  imparted  velocity  in  Zg  direction 
P*  = pressure  at  muzzle  for  critical  conditions 
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+ n = orientation  of  asymmetric  force  in  the  muzzle 
jet 

= plan  area  of  one  fin 
Mp  * mass  of  projectile 
The  transverse  angular  velocity  is 


where 


i' 


3 + i 

o 


/ _ 


= (Y+l)  p 2A 


AD  Z 

f 2 

1 I V i 
Y P 


P (e 


i C4>  + 

e VTc 


)>  (2) 


I = transverse  moment  of  inertia 

y 


A * moment  arm  of  force  on  fins  about  the  center 
of  gravity. 

Finally,  the  reverse  spin  imparted  by  the  muzzle  blast  due  to  differ- 
ential fin  cant  with  an  equivalent  moment  arm  length  rQ  is 


x 

♦' ' = - (Y+l)  p n Af  ° D \ p 6 (3) 

0 1 I V Z 

x p 

where 

I = axial  moment  of  inertial 
x 


6 = angle  of  differential  fin  cant 


III.  THE  JUMP  EQUATIONS 

The  dynamic  state  of  the  projectile  upon  leaving  the  muzzle  blast 
region  and  entering  into  free  flight  defines  its  subsequent  trajectory 
permitting  the  computation  of  jump  due  to  muzzle-gas  effects.  The 
deviation  from  its  intended  trajectory  is 

Y + i Z 

Q = lim  — (4) 

Y 

X -*•  co  e 
e 
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where  the  coordinate  system  is  given  in  Figure  3.  Here  the  influence 
of  gravity  is  not  considered. 

Murphy  and  Bradley1  show  that  the  aerodynamic  jump  is  approximately 


°j  - Jf'  f'o  * JA  <VO  CS) 

The  first  term  on  the  right  hand  side  is  the  jump  due  to  the  initial 

transverse  angular  momentum  given  the  projectile  by  the  muzzle  blast. 

For  the  asymmetric  projectile,  £ / is  given  by  Equation  (2).  The 

0 1 

expression  for  the  coefficient  in  the  first  term  is 


V 


(6) 


M £ 
P 


“M 


While  the  first  term  in  Equation  (5)  is  applicable  to  both  symmetric 
and  asymmetric  projectiles  launched  with  an  initial  angular  velocity. 


the  second  term  involves  only  the  asymmetry  of  the  projectile.  Here 

C.  » C.,  are  the  free  flight  lift  and  static  moment  coefficients, 
a 

The  expression  for  is 


pS£ 
A 2 M 


) e e 


<t> 


(7) 


and 


where 


> 


C = 


initial  free  flight  roll  rate.  Equation  (3) 
steady  state  roll  rate* 

- \ {/ictp  * cyy2>  v <8> 

roll  damping  coefficient1 

- [pS£/2M  ] t(Ix/M £2)  + CD]  (9) 

P 

roll  moment  coefficients  due  to  roll 
and  fin  cant  respectively 
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The  expression  for  will  be  given  and  discussed  in  detail  in  the  next 
section. 

To  calculate  the  total  deviation  of  the  projectile  from  its  intended 
trajectory,  one  must  also  consider  the  transverse  linear  velocity  im- 
parted to  the  projectile  by  the  muzzle  blast.  Adding  this  effect  to  the 
aerodynamic  jump,  one  obtains  for  the  total  deflection: 


0 * vvP  + jc7  + (JA  V*'  -)  • 

Using  Equation  (1)  and  (2),  one  obtains 


(10) 


0 ■ [1  + j ] (yl)  p*  2A  

\ MpVP 


= , i ( <()  +Tr), 
P (ee  VYe  J] 


+ W O • 


(11) 


IV.  APPROXIMATION  OF  * 


The  expression  for  $ is 


t'  st  % 

*t  = lim  ^Cs-'1  0 fo  exp  ^ ® dst  dst  (i2) 

_ . „ t 12 

st 

where  the  roll  angle  is  given  as 

♦ “ (to'/C)  (st-  [ / 4 »„')  - 1]  texp(-st)-l]}  , (13) 

and 

st  = 0S 

s = distance  along  the  trajectory  in  calibers. 

In  this  section,  a series  approximation  to  $t  (A,B)  is  given  in 
terms  of  the  parameters: 

A 5 ij/c  , 

» = ■ 
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Murphy  and  Bradley1  have  shown  that  the  integral  in  Equation  (12)  may  be 
transformed  to 

* (A,B)  = iA  ! ™ exp  [Af(r)]  dr,  (14) 

where 

f (r)  = - (r  + i (B— 1)  [exp  (-ir)-l]}  . (15) 

They  further  obtained  an  expression  for  * (A,0)  when  A is  large: 

(A,  0)  = [(1+i)  (irA)^/2)  + (i/3)  (16) 

For  the  current  application,  it  is  necessary  to  consider  small, 
but  non-zero,  values  of  B;  thus,  <J>^  (A,B)  is  expanded  in  a Maclaurin 
series : 

d4>  d2$ 

4>r(A,B)  w *t(A,0)  + -*g-  | B + h | B2  + ..  (17) 

B-0  dBZ 

B=0 

Differentiating  4>  with  respect  to  B,  we  obtain 
d4>  » 

= A""  [exp(-ir)-l]  exp  [Af(r)]  dr  (18) 

Differentiating  4 again,  we  obtain 

d2$t  3 “ 2 

— = - i A * f [exp(-ir)-l]  exp  [Af(r)]  dr  (19) 

dB  o 

For  large  A and  small  B,  the  main  contribution  to  the  integral  in 
Equation  (14)  will  occur  when  r is  small.  When  r is  large,  the  integrand 
consists  of  small  oscillations  about  zero,  which  tend  to  cancel  out  in 
the  integration.  For  B=0,  Equation  (15)  may  be  expanded  in  series  form: 


f (r)  = - (ir2/2)  - (r3/6)  + (ir4/24)  + (r5/120)  + ...  (20) 
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Using  this  series,  we  can  write 


exp  [Af(r)]  ■ [exp(-iAr2/2)]  exp  (Au),  (21) 

where 

u = - (r3/6)  + (ir4/24)  + (rS/120)  + (22) 


In  turn,  the  term  exp  (Au)  is  expanded  in  series.  While  Murphy  and 

Bradley  needed  to  consider  only  two  terms  in  the  expansion,  Equation  (21), 

to  obtain  their  approximation  for  4>  (A,0),  the  present  authors  carried 

9 1 

the  expansion  out  to  the  r term  for  their  approximation  to  ^(AjB). 
Finally,  expanding 


exp(-ir)  - 1 *»  - ir  - (r2/2)  + (ir3/6)  + r4/24,  (23) 

Equations  (21)  and  (23)  may  be  substituted  into  Equations  (17),  (18), 
and  (19).  By  a transformation  of  variable,  we  obtain  integrals  of  the 
form  „ 

oo 

/ rim  exp  (-ip2)  dn 

o 

Integrating  this  integral  around  the  contour  which  includes  the 
positive  part  of  the  real  axis  and  the  line  45°  from  the  real  axis, 
we  obtain 

/ n"  exp(-in2)  dn  » exP  r (2lL) 

0 

where  r is  the  well-known  gamma  function. 

By  evaluating  the  gamma  function  and  neglecting  negative  powers 
of  A,  * is 

*t  = ( [(1+i)  (-rrA)Js/2]  + [i/3]}  + [-A  + (2i/3)]B 

+ { [ (1-i)  (irA)^  A/4]  - [A/3]  + [ (1+i)  firAD^/48]  + [2i/135])B2.  (24) 

This  approximate  expression  for  4>  is  compared  with  the  exact  numerical 
2 r 

computations  in  Figures  5 and  6.  The  approximation  is  quite  good  up 
to  B = ± 0.1  but  fails  for  greater  absolute  values  of  B as  A becomes 
large. 
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Figure  5.  Magnitude  of  4>t  versus  <{£  C-"* 
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100 


Figure  6.  Argument  of  4>t  versus  C 


In  Equation  (24) , it  is  apparent  that  the  approximation  becomes  worse 
for  large  A since  the  leading  term  in  the  factor  that  multiplies 
Bn  is  The  limitation  does  not  affect  the  analysis  since  B 

will  seldom  have  a magnitude  greater  than  0.1.  For  instance,  the 
value  of  B found  for  a 5.8mm  flechette  launched  from  a smoothbore  tube 
is  only  -0.0403. 

V.  DEVIATION  CAUSED  BY  MUZZLE  BLAST  ON  ASYMMETRIC  PROJECTILES 

The  total  deviation  of  the  projectile  is  given  by  Equation  (11); 
however,  even  if  the  projectile  could  be  launched  in  an  environment 
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free  of  muzzle  blast,  it  would  still  experience  a jump  due  to  asymmetry 
equal  to  (A,0)/ ] . Thus,  the  contribution  of  muzzle  blast  to 

the  total  deviation  is  defined  to  be: 


°mb  = 0 - IJA  «t  ’ (25) 

From  Equation  (24),  the  lead  terms  of  4>^  - $ (A,0)  are 


*t  - *t  (A,  0)  * - AB  { 1 — [ 2i/  (3A)]  - [ (1-i)  (uA)**  B/4] } 

= - ABT.  (26) 

Thus,  from  Equation  (7)  and  (11) 

C, 


®mb  '[“Cl1  ■>*  2 V ^T2 


P 


M 


PS1  c ,,  Cl»  \ 

M ST  CN 

a e 


2M  N 
p e 


M V 
P P 


x i 4>  ABT 

) ee  e 

~cr. 


(27) 


According  to  linear  theory, 

C 


M A 

E _ C 

:N  =‘*  ’ 


where 


= c.g.  - c.p.  separation  of  asymmetric  fin 
~ A. 


Substituting  into  Equation  (27), 


0mb  = - [1  + p— - £]  ee1^  {(y+1)  p 2A 
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CM  * 
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f 2 

MV 
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. PS i r ABT  , 
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(28) 
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Examine  the  right  hand  term  within  the  braces  in  Equation  (28): 


- pSt  r 

~ 2M  N 
P e 


ABT 

<j/ 

T00 


psi  r l \ 

2M  N C j 

p E c 


(29) 


From  Equations  (8)  and  (9), 


C<$f  = C.  6 
" *6 

yielding: 


M 6Jt2 
P 

Again,  according  to  linear  theory 


Substitution  of  Equations  (3)  and 


2 (pS£)/2Ix), 


2 

2 l . 
n r 

o 

(32)  into  Equation  (31)  yields 


(30) 


(31) 


(32) 


¥ - ~ Cy+1)  P 2 Af  --p-  a P T.  (33) 

1 M V Z 
P P 


Equation  (28)  becomes: 

CL 

®mb  - - < t1  * cf  I’  (Y*1)  P’  2Af  77^  P 1 [1-T>  Cel  (34) 

M MV 

a P P 

The  quantity  in  the  left  hand  brace  is  identical  to  the  coefficient 

2 

in  the  deflection  equation  describing  the  muzzle-blast -induced 
deviation  of  a symmetric  projectile  launched  at  an  angle  of  attack,  5 . 

TT1 

Hence,  by  evaluating  the  magnitude  of  1-T,  a direct  comparison  can  be 
made  between  the  sensitivity  of  symmetric  and  asymmetric  projectiles 
to  muzzle  blast.  From  Equation  (26), 

1 - T = [2i/  (3A)  ] + [(1-i)  (ttA)5*  B/4)  . 
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The  value  of  A is  frequently  near  25  while  B is  on  the  order  of  0.05; 
thus,  the  magnitude  of  1-T  is  approximately  0.1.  Since  a typical 
muzzle  blast  deviation  of  a symmetric  projectile  is  on  the  order  of 
0.1  mil  deflection  per  degree  of  initial  angle  of  attack,  the 
sensitivity  of  the  asymmetrical  projectile  would  be  about  0.01  mil 
deflection  per  degree  of  asymmetric  fin  deflection. 

VI.  SUMMARY  AND  CONCLUSIONS 

Using  a simple  approximation,  a closed  form  expression  for  the 
effect  of  muzzle  blast  on  the  jump  due  to  asymmetry  of  fin-stabilized 
projectiles  is  obtained.  The  analysis  shows  that  muzzle  blast  does 
not  significantly  alter  this  portion  of  the  total  jump. 
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LIST  OF  SYMBOLS 


planform  area  of  single  fin 

local  sound  speed 

roll  damping  coefficient 

roll  moment  coefficient  due  to  roll 

roll  moment  coefficient  due  to  fin  cant 

projectile  lift  coefficient  and  derivative  with  respect 
to  a 

projectile  static  moment  coefficient 

moment  coefficient  due  to  fin  deflection 

normal  force  coefficient  due  to  fin  deflection 

distance  from  obturator  to  fin  c.p. 
diameter  of  gun  bore 

longitudinal  and  transverse  moments  of  inertia 

jump  coefficients 

shaft  diameter  of  projectile 
dimensionless  fin  lift  function 
Mach  number  or  moment 
mass  of  the  projectile 

number  of  fins 
pressure 

dimensionless  momentum  function 
muzzle  jet  and  in-bore  momentum  impulse 

distance  along  trajectory  in  calibers  (s  * x /£) 

6 

2 

reference  area  of  projectile  (S  = ir£  /4) 

time  (t=0  corresponds  to  obturator  passing  muzzle) 


w 

X,Y,Z 

A#  ^ 

a,  6 

Y 

6 

A,  Af 


e 

£ 

0 

©. 

1 


mb 

P 


Subscripts 


e 

o 


Superscripts 

o 
( )' 

( )* 


LIST  OF  SYMBOLS  (Continued) 
projectile  velocity 

transverse  velocity  in  Zg  direction 
coordinates 

angles  of  attack  and  sideslip  in  non-rolling  coordinate 
system 

ratio  of  specific  heats 
differential  fin  cant  angle 

c.p.  - c.g.  separation  in  reverse  and  forward  flow, 
respectively 

magnitude  of  asymmetric  fin  cant  angle 

At  A# 

complex  angle  of  yaw  $ + i a 

angular  deflection  of  projectile  from  boreline 
aerodynamic  jump:  angular  deflection  from  particle 

trajectory  (gravity  and  drag  determined)  due  to 
aerodynamic  forces 

angular  deflection  of  projectile  due  to  muzzle  blast 

density 
roll  angle 

initial,  free  flight  orientation  of  asymmetric  fin  force 


denotes  earth-fixed  coordinates 
denotes  projectile  properties  immediately  after 
penetration  of  the  muzzle  blast 
denotes  ambient  or  steady  state  conditions 

denotes  dimensionless  quantities 
denotes  derivative  with  respect  to  s 
denotes  critical  or  sonic  conditions 
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